Neutrinos decoupled from /3-processes and supernova explosion 
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Based on the gravitational collapse time-scale is larger than the weak interaction time- 
■ scale at core densities p > 10 n gr/cm 3 , we approximately use the /3-equilibrium condition and 

(N 

. ■ particle number conservations to calculate the number and energy densities of neutrino sphere 

Ch ■ 

in the process of gravitational core collapse towards the formation of a proto-neutron star. We 
, find that at core densities pjee > 10 12 gr/cm 3 , the /^-equilibrium condition cannot be satisfied 

consistently with charge, baryon and lepton number conservations, leading to the presence 
of excess neutrinos decoupling from the /3-equilibrium. These excess neutrinos interact with 
nucleons and electrons via the neutral current channel only and their diffusion time is about 
10 2 sec. The excess neutrino flux could play an important role in an Supernova explosion, 
O , provided the fraction of excess neutrinos over all neutrinos is at least one present, 

to 

PACS numbers: 13.15.+g, 97.60. Jd, 26.30. Jk 
Introduction. 

A great effort has been made to understand the phenomenon of a core-collapse supernova for 
a long time. It is known that the dominant weak interaction process altering the composition of 
q \ the core matter in during collapse is /3-process (electron capture and neutrons decay) on free and 

bound nucleons, which proceeds at a rate sufficient to produce a large number of neutrinos. In 
the core at densities larger than 10 11 gr/cm 3 , these neutrinos are trapped and thermalized, leading 



X 



to electron capture^ equilibrium over time scales shorter than the characteristic dynamical time 

The collapse continues essentially homologously [4], until nuclear densities 



1U 



scales for collapse 

reach 10 14 gr/cm 3 . Since nuclear matter has a much lower compressibility, the homologous core 
decelerates and bounces in response to the increased nuclear matter pressure. This drives a shock 
wave into the outer core, i.e. the region of the iron core which lies outside of the homologous core 
and in the meantime has continued to fall inwards at supersonic speed . If the shock wave were 
to propagate outward without stalling and make an explosion with energy about ~ 10 51 ergs, it 
would be a successful prompt hydrodynamical explosion, but all of the realistic models completed to 
date suggest that this mechanism does not occur at least for massive collapsing iron core, because 
the shock wave loses energy in dissociating iron nuclei when this shock passes through outside 
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matter (the outside core matter includes iron group nuclei). Then the shock wave is enervated and 
loses its energy in the form of electron neutrinos, finally the shock wave would halt its outward 
motion [IflE]. 

Wilson proposed the delayed mechanism || that neutrino flux has a time scale much longer 
than prompt hydrodynamical explosion and revives the stalled shock wave by the charged current 
absorption of electron neutrinos and anti-neutrinos. However, in this mechanism, at least 99% of 
the binding energy of the neutron star (~ 10 53 ergs) comes out in neutrinos, which is 100 times that 

n a 

needed for the shock wave to give a powerful supernova explosion [5J, ]6|] . In addition, compared 
with the time scale of the shock wave, the diffusion time of these neutrinos is too much longer 
to revive the shock wave at proper time jjjj]. The main problem of delayed mechanism is then 
channeling some small fraction of the neutrino energy to the proper place and at the proper time 
to cause the explosion, several attempts have been made to find a solution to the problem [g, |7J, W\ ■ 
In this letter, we study a possible solution to the problem. 

Weak interaction and adiabatic gravitational collapse. In the standard model of particle 
physics, neutrinos interact with electrons and nucleons via charged and neutral current processes, 



e +p =f± n + u e , 
z/+(e, N) -»■ v + (e, N). 



(1) 
(2) 



The cross-section of the dominate /3-processes ([I]) a cc ~ 9.75 x 10 42 (Ev/10MeV) 2 cm 2 where E u 
is the neutrino mean-energy. Neutrinos are left-handed, interacting only with left-handed quarks 



(u, d) inside nucleons and electrons by exchanging charged gauge bosons 



V^. The cross-section of 



the neutrino-electron channel ([2]), a 
the neutrino-nucleon channel ([2]) is 



ei> 



0.01o- cc for E u ~ 10 MeV 
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111 ] . The cross-section of 
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(3) 



where g u L = 1/2 - (2/3) sin 2 W , g d L = -1/2 + (1/3) sin 2 W , g u R = -(2/3) sin 2 W and g d R = 
(1/3) sin 2 9\v are left-handed neutrino gauge couplings to left- and right-handed quarks (u,d) by 
exchanging neutral gauge boson Z° (sin 2 9w ~ 0.25) U]. <t nc w a cc /A, the /3-processes ([I]) have 
a larger probability than neutral current processes ((2j). 

We consider a collapsing stellar core of radius R and mass M ~ M@, the collapsing time-scale 
icon = (R/R) ~ R/ c is much larger than the weak-interaction time scale i wca k ~ ( ccr cc n «) _1 
(*coli > iwcak), at nucleon densities p > 10 11 gr/cm 3 for R < R* sa 20R C and E v ~ 10 MeV (see 
Fig. [1]). This implies that gravitational collapse could be approximately treated as slowing varying 
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FIG. 1: In the range -Rdcc < R < R*, the ratio i we ak/icoii is plotted as a function of the collapsing radius 
R/Rc- 

adiabatic process with respect to the /3-processes ([I]), which is in agreement with other studies (see 
for review jl^]). In this adiabatic approximation, at each collapsing radius R, the /3-equilibrium 
is assumed to be locally and instantaneously established in the characteristic space-time variations 
AR and At, determined by the equation of gravitational collapse. The analogous discussion can be 
applied for other macroscopic processes. Based on this adiabatic approximation, we study neutrino 
emission in gravitational collapse. 

/3-equilibrium and neutrino sphere. Suppose that the stellar core is composed by complete 
degenerate gases of electrons, neutrinos, protons and neutrons for their Fermi momenta being much 
larger than temperature. The /3-equilibrium condition is 

Hn + V-v = Hp + He, (4) 

where chemical potentials 



Hn,p,e,u — \J ( c Pn,p,e,u) 2 + m n,p,e,u ci ■ (5) 

Suppose that particles homogeneously distribute within the core of volume V = 47ri? 3 /3, their 
densities n n>p>eil/ = N n ^ e<v /V and Fermi momenta Pn, p , e ,u = i^ 2 ) 1 ^ n n,p,e,w A = N n + N p 
and L = N e + N u are the total conserved baryon and lepton numbers, local neutrality requires 
n p = n e . The variation of neutrino number is related to nucleon one, AN U = AN n = —AN p . These 
equations completely determine the ratio A/N p and n n)P)( , yU (R) as functions of the core radius R. 
During collapse, the variation of gravitational binding energy A£ g ~ — 3(GM 2 / R) (AR/R) [jij, 



and it is believed that about 99% of this energy is converted to the neutrino energy 



141 ] . Neutrinos 
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FIG. 2: In the range p* < p < pdec: the ratio (A/Np) 1 / 3 is plotted (solid line) as a function of the core 
density p. The dashed line is for the case all neutrinos escape p v = 0. 

are trapped inside the core for large opacity 

t v ~ {(J cc +a* c )n n R&a 0C n n R> 1, (6) 

for R < R* ~ 20R C , coinciding with the valid range of adiabatic approximation. 

Using Eqs. (JU[6]), we determine a neutrino sphere (see Figs. [2] and [3]) by numerically calculating 
the ratio A/N p , particle number-density n n ^ v ^ e ^ v and energy-densities e n ,p,e,v as functions of the core 
radius (or core density), starting from the radius R* ~ 20i? c (p* ~ 5 x 10 10 gr/cm 3 ) to the radius 
Rdec ~ (2 ~ 5)i? c (pdcc ~ Pnuci(^c/-Rdoc) 3 )i which will be clarified below. The obtained value 
A/Np at the radius -Rdec ~ Rc ( a new born proto- neutron star radius) is consistent with numerical 
simulation results Y" e ,p = ^e,p/A « 0.3 — 0.5 0, [3j in agreement with observation data {l^ ]. 

It is shown that at the radius -Rdeo T ^ ~ 10 2 ~ 3 , n v ~ 10 36 ~ 37 /cm 3 , ej/ ~ 10 31_32 ergs/cm 3 , and 
the total neutrino energy and number are 10 52_53 egrs and 10 56-57 . The neutrino mean-energy 
E u = e u jn v ta lOMeV. 

Excess neutrinos over the /3-equilibrium. The numerical calculation shows that at the radius 
-^dec ~ (2 ~ 5)i? c , where the neutrino chemical potential is so large that non solution fully fills the 
/3-equilibrium condition @ with baryon, lepton and charge conservations. In fact, it is a necessary 
condition that the neutrino chemical potential should be smaller than electron one, \i v < /i e , 
because the neutron chemical potential is larger than proton one fj, n > fi p . In another word, 
N v < N e /2, where the factor 1/2 is due to different neutrino and electron spin degeneracies. We 
might introduce the radius i?d ec determined by the critical condition fjL u (Rdec) = ^e(Rdec), then 
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FIG. 3: In the range p* < p < p^cc and p& cc re p n uci(Rc/ Rdec) 3 where i?dec — (2 ~ 5)i? c , we plot the 
homogeneous number-density n v (left panel) energy-density t v (right panel) of neutrino sphere formed at 
the core density p. 



Eq. (j4]) yields the critical condition 

N n ~ 2iV„ = N e = N p = A/2, (7) 

at the radius i?dcc- We analytically calculate the neutrino chemical potential fJ,v(Rdec) = 
/j, e (Rdec) ~ 250MeV(i? c /i?dcc), number-density n l/ (Rd cc ) — 2 x 10 38 (i? c /i?dec) 3 cm -3 and energy- 
density e v (i?dcc) — 2 x 10 34 (i? c /i?dec) erg cm -3 , consistently with the numerical result at the radius 
i?dcc (see Figs. [2] and [3]). 

When the collapsing radius R < Rdec: N v > N n /2 then neutrons cannot absorb all neutrinos 
via the /3-process n + v —> p + e. As a consequence, excess neutrinos N„ X over the /3-equilibrium, 
namely those neutrinos decoupled from the /3-processes (JT]), must be present in the neutrino sphere. 
This can be understood from the microscopic point of view, the difference between the neutrino 
emissivity and the absorption by the thermal system of electrons, protons and neutrons via the 
/3-processes ([TJ is given by [ijj 

Cp = J du v [B v - J„), Qp = J du u -^y (By - J„), (8) 

where k u , I u and B v are the absorptive opacity, the specific intensity and the black-body function 
for neutrinos. In Eq. ©, J 7 ^ = [ e K--Oe+^ P -Mn)]/feT + anc j w ^ ^ | p ^|_ "when the detailed 
balance is established, i.e., the neutrino absorption and emission rates are exactly equal, neutrinos 
are in thermal equilibrium with the system via the /3-processes, the /3-condition ([I]) is fully satisfied 
and C/3 (Qp) vanishes, leading to the neutrino black body distribution T v = T v = [e( Wl/ ~^ l/ )/ feT +l] _1 
and By = I u . This is the case for R > i?dec , we use /3-equilibrium condition dU) to calculate number- 



0.6 r 



0.5 - 



0.4 - 



0.3 - 



0.2 - 



0.1 



0.0 L 



lOp* 



lOOp* 



~10 3 p* 



FIG. 4: Y v (p)/Y e (p) is plotted as function of core density (* and quadrangle points are results of [18[ and 



fitting formula 



19(), p* ~ 5 x 10 10 gr/cm 3 is the density which is started the trapping of neutrinos. 



and energy-densities of the neutrino sphere. In the case R < i?deo the /3-equilibrium condition 
([T]) cannot be satisfied and T v ^ J-' u , and the neutrino chemical potential /u£ on is different from 
the /3-equilibrium one /i„ (piJJ on > fi u ). As a result, Cp, Qp / indicates neutrinos decoupled from 
the /3-processes ([1]), namely the excess of neutrinos over the /3-equilibrium. We call the radius 
R&ec ~ (2 ~ 5)i? c decoupling radius, corresponding core densities pdec ~ (10Km/i? dec ) 3 p nuc i, where 
Pnuci is the nuclear density and R dcc = (20 ~ 50)Km. 

We would like to see whether this critical condition ([7]) appears in calculations by numerical 



simulations. Using numerical simulation data 



181 ] and fitting function [19(] for electron fraction 



Y e , in Fig. @, we plot the ratio Y u /Y e (neutrino fraction over electron fraction) in terms of core 
density p. We find at core densities ~ 10 3 this ratio Y u /Y e — >• 0.5, consistently with the decoupling 
density pdcc obtained in Fig. [2] and Eq. ([7]). 

Excess neutrino luminosity and diffusion time. These excess neutrinos decouples from the 

EX 

/3-processes, then the cross-section of excess neutrinos via the /3-processes vanishes, a cc = 0. This 
implies that all left-handed quark fields inside nucleons are fully engaged with left-handed neutrinos 
via the /3-processes when R < -Rdcc- Therefore, excess neutrinos scatter only with right-handed 
quarks via neutral current interaction, i.e., the g^' d -terms in Eq. ([3]), yielding the cross-section 
°nc ~ 0.01<7 CC . As a result, the excess neutrino opacity is only due to the scattering via right- 
handed neutral current channels ([2]), 



EX EX t> EX ... T> 

t v = V *;R~0Mt u , * m =o% 1 +o* (9) 
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where ~ n N /2 the nucleoli number-density, which accounts only right-handed quark fields inside 



nucleons. The diffusion time t^ s can be estimated by (see for example 12]), 



C - 3^iV coll ~ 3^£ ( Rg) = 3 ^££ T f, (10) 



EX 



where the collision number iV co n ~ {Rdec/^u ) > an d the mean-free path of excess neutrinos Aj, 
l/(o"^n^). At the radius Rdea we obtain the excess neutrino diffusion time t^* s ~ 0.01 sec and 
averaged diffusion velocity v^q = Rdec/t^ts ~ 10~ 1 c. The outgoing energy-flux of these excess 
neutrinos F v = Qp = r/ EX e u v^ s , where the parameter ?? BX (r? EX < 1) presents the fraction of excess 
neutrinos over the /3-equilibrium. The corresponding luminosity is 

r EX 2 EX 

L u (R) = 4ttR t? ex e v v dm 

( \ EX \ / R \ 2 



= 5 x 10 57 r/ EX ^ j j erg/sec, (11) 
and L l/ (i?dec) ~ 5 x 10 54 r/ EX erg/sec, increasing as the collapsing radius R decreases in the range 

Rc < R < Rdec- 

We try to estimate the fraction ry EX of excess neutrinos produced in the range R c < R < i?dec- 
Due to their short diffusion time ~ 10 -2 sec, excess neutrinos are supposed to completely emit 
away when the collapsing radius R varies from Rdec to the radius R c of a new born proto-neutron 
star, in which protons, neutrons, electrons and /3-neutrinos are in the /3-equilibrium at the critical 
condition N n ~ N p = N e = 2N U . Under this assumption, we calculate the energy variations 

AEp%,e,v = E p%,e,v(Rc) - E™n,eA Rdec ^ — (°- 9 ~ x l° 53er g s i as wen as the energy emission 
of excess neutrinos (jlip . The variation of gravitational binding energy A£ g ~ —3(GM 2 / R)(l — 
Rc/ Rdec) ~ — (1 ~ 1-7) x 10 53 ergs. We obtain rj EX < 0.1 by the total energy conservation. 
Excess neutrinos and Supernova explosion. The necessary conditions for an Supernova 
explosion due to neutrino flux are: (i) a steady neutrino flux occurs in the range between a proto- 
neutron star (PNS) radius R c (shock wave starts) and shock radius -R s hock ~ 200km (shock wave 
stalls), corresponding to the shock time-interval i s h ~ 0.1 sec; (iii) the steady neutrino luminosity 
exceeds a threshold luminosity 1 . The theoretical and numerical studies show in order to have a 
powerful explosion, it needs the neutrino luminosity threshold L" A ~ 5 x 10 52 (M/Mq) erg/sec to 
overcome the gravitational pressure of infilling matter. The rate of infilling matter M/Mq ~ 0.1 
decreases during the shock time interval, and increases after shock wave stalls [17I ]. If the neutrino 
diffusion time is smaller than or the same order as the shock wave time, the outgoing neutrino flux 
is in the same direction of outgoing shock wave with small rate M/Mq. This implies that (i) the 
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critical luminosity L"" becomes small; (ii) neutrinos have much more probability to interact with 
infilling matter and transfer their energy-momenta to make a powerful explosion [71. 191. 1171] . 




We are in position of marking the following remarks on the excess neutrino luminosity (jlip for 
an supernova explosion: 

1. The excess neutrino luminosity turns on at the radius Rd cc ~ (2 — 5)R C , indicating that an 
instability is triggered. The shock wave occurs around this radius; 

2. The time scale of excess neutrino luminosity is about 10~ 2 sec, the same order of shock 
wave time scale t s h, indicating that the excess neutrino luminosity starts to act for explosion 
almost at the same time as the shock wave starts, rather than after shock wave stalls; 

3. If the fraction of excess neutrinos overall neutrinos is about one percent (r/ EX ~ 1%), the 
excess neutrino luminosity ([lip is larger than the threshold luminosity L"^ ~ 5 x 10 52 erg/sec, 
satisfying the necessary condition for an supernova explosion. 

These properties of excess neutrino luminosity might give a solution to the main problem of delayed 
mechanism: channeling some small fraction (r/ EX ~ 10~ 2 ) of the neutrino energy to the proper place 
[Rdcc ~ (2 — 5) x 10 6 cm] and at the proper time (i^ iff ~ 10 -2 sec) to cause the supernova explosion. 
The energy of excess neutrinos can be estimated by L U {R& CC ) x ~ ?/ EX 10 53 ergs. However, it 
should be mentioned that (i) the excess neutrino decoupling radius .Rdec ~ (2 — 5)R C is a low- 
limit, because there are other processes of neutrino productions which we do not consider; (ii) 
the excess neutrino diffusion time-scale 10 -2 sec is also a low-limit because it is obtained from 
the purely right-handed neutral current interaction; (iii) the excess neutrino fraction r/ EX < 0.1 is 
an up-limit, because it is obtained by total energy conservation and we do consider all possible 
energy dispassive channels. In the present model and calculations, we are not able to give the up- 
limit of the excess neutrino diffusion time-scale and the low-limit of excess neutrino fraction r/ EX . 
Thus, excess neutrinos might not play an important role in Supernova explosions for the following 
two cases: either the the excess neutrino diffusion time-scale is larger than the shock time-scale 
(tth ~ 0.1 sec) or the the excess neutrino fraction is so small that the excess neutrino luminosity is 
smaller than critical one L?."*. Nevertheless, we will show that in cooperation with strong electric 
fields on stellar core surface [201 ]. these excess neutrinos can at least play an important role for 



Neutrino emission in cooling phase. In the neutrino sphere, in addition to excess neutrinos, 
most neutrinos participate the /3-equilibrium ([I]) and we call them /3-neutrinos to distinguish them 




enhancing electric field 




and trigging electron-positron pair productions [221 ]. 
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from excess neutrinos. Their diffusion time can be estimated by 



£ diff - 3^prJ, tI = n N (a cc + <r NC ) « n N a cc . (12) 



which is about a few seconds (t^ iff > lsec) for the core radius R c ~ 10 6 cm and opacity t v ~ 10 3 . 
The diffusion time ()12|) of /3-neutrino flux is much larger than the shock time t s h, the outgoing 
/3-neutrino flux starts after the shock wave stalls and the infilling mater rate M/M & increases. 
This implies that (i) the threshold luminosity L" 1 becomes large; (ii) /3-neutrinos have not much 
probability to interact with infilling matter and transfer their energy-momenta to make a powerful 
explosion. Therefore /3-neutrinos cannot be relevant for the delayed mechanism to revive the shock 
wave. 

Because of their short diffusion time (10 _2 sec), excess neutrinos are assumed to completely 
emit when a new born proto-neutron star is formed at the radius R c after an Supernova explosion. 
The number of major /3-neutrinos left over inside the neutrino sphere is about A/2, which is 
approximately given by the critical condition ([7]). Due to the radial gradient of pressure and 
particle density inside the core, these /3-neutrinos diffuse outwards and emit from the core, leading 
to the cooling of new born proto-neutron stars. In this cooling phase, while /3-neutrinos are diffusing 
out of the core, the neutrino production via the process of electron capture ([I]) takes place. As 
a consequence, the value A/N p increases from the value A/N p rs 2 for neutrinos trapping to the 
value A/Np ~ 250 for all neutrinos escaping at the end of cooling phase (see Fig. [2]). Thus, the 
total number of neutrinos produced by the process of electron capture should be about A/2. The 
total number of emitted neutrinos in the cooling phase should be about A, which is the sum over 
/3-neutrinos and neutrinos produced by electron capture process from A/N p ~ 2 to A/N p ~ 250. 
Carrying major neutrino energy and number in the entire neutrino sphere, /3-neutrinos play an 
important role in the cooling phase of a new born proto-neutron star ja,[3|. As mentioned, excess 
neutrino energy is about r/ EX 10 53 ergs. Then /3-neutrino energy should be about (1 — ?7 EX )10 53 ergs. 
This is the predominant cooling mechanism immediately after formation proto-neutron star with a 
timescale of seconds, and the total energy carried away by /3-neutrinos and neutrinos produced by 
electron capture process should be the major part (about 99%) of the gravitational binding energy 



B,Q. 



Some remarks. In this Letter, we only consider the core of collapsing star at radii R < R* ~ 
20R C (p* ~ 10 11 gr/cm 3 ) on the basis of homologous collapse [4]. However, we assume that the 
core density p is spatially homogeneous and its sharp boundary described by a Heaviside step 
function, i.e., p = 5{r — R)M/V , neglecting the variation of core density at the boundary. This 
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approximation leads to the effective mass inside the core is larger than realistic one, then energy- 
and number-densities of neutrino sphere calculated are probably larger than that in realistic case. 
In addition, ignoring nuclei of iron group in the core, we adopt a model of completely degenerate free 
electrons, protons and neutrons to calculate their chemical potentials. This approximation affects 
on the calculations of the ratio A/N p and neutrino productions. Under these approximations, we 
show the possible presence of excess neutrinos decoupled from the /3-processes, and obtain the 
preliminary results of the decoupling radius Rdea excess neutrino fraction 7? EX , excess neutrino 
diffusion time t^jj and excess neutrino luminosity , which could be a possible solution to the 
main problems of delay mechanism for supernova explosion. We also briefly discuss /3-neutrinos 
and neutrinos produced by electron capture process in the cooling phase of new born proto-neutron 
stars. The detailed calculations to obtain these preliminary results of this Letter will be presented 
in a lengthy article 



221 ] . Needless to say, these preliminary results are necessary to be verified by 
other approaches, in particular numerical approaches. 
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